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Abstract 

We study the factorial divergences of Euclidean a problem with connec¬ 
tions both to high-energy multiparticle scattering in d=4 and to d=3 (or 
high-temperature) gauge theory, which like is infrared-unstable and su- 
perrenormalizable. At large external momentum p (or small mass M) and 
large order N one might expect perturbative bare skeleton graphs to behave 
roughly like N\{ag‘^/p)^ with a > 0, so that no matter how large p is there 
is an A ~ pjg^ giving rise to strong perturbative amplitudes. The semi- 
classical Lipatov technique (which only works in the presence of a mass) is 
blind to this momentum dependence, so we proceed by direct summation of 
bare skeleton graphs. We hnd that the various limits of large A, large p, 
and small mass M do not commute, and that when A S> p‘^/M‘^ there is a 
Borel singularity associated with g'^/M, not g'^/p. This is described by the 
zero-momentum Lipatov technique, and we hnd the necessary soliton for 
the corresponding sphaleron-like solution for unbroken Yang-Mills (YM) the¬ 
ory has long been known. We also show that the massless theories have no 
classical solitons. We discuss non-perturbative effects based partly on known 
physical arguments concerning the cancellation by solitons of imaginary parts 
due to the perturbative Borel singularity, and partly on the dressing of bare 
skeleton graphs by dressed propagators showing non-perturbative mass gen¬ 
eration, as happens in d=3 gauge theory. 
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1 Introduction 


Several interesting questions arise about the large-order behavior of pertur¬ 
bation theory in (Euclidean) infrared-unstable theories like d=3 QCD or elec- 
troweak (EW) SU{2), which have dimensionful coupling constants and are 
superrenormalizable. These theories have physical interest in themselves as 
the iV = 0 Matsubara-mode magnetic sector of high-temperature (above the 
phase transition) gauge theories, which is strongly-coupled even for EW the¬ 
ory, and there is a connection also to the behavior of large-multiplicity high- 
energy scattering amplitudes. The Wilson-Fisher epsilon-expansion does not 
work for the three-dimensional gauge theories (there is no hxed point), and 
in these theories, massless in perturbation theory, there is no classical sad- 
dlepoint of hnite action to use with the Lipatov[|I| technique, the standard 
method for dealing with large-order perturbation theory. Even if there were 
such a classical saddlepoint this well-known technique would fail in the semi- 
classical approximation, and the failure is precisely of the type encountered in 
the Ringwald-Espinosa0 attempt to hnd strong baryon violation in the EW 
sector]^]. Ringwald and Espinosa were concerned with on-shell high-energy 
large-A^ amplitudes in Minkowski space|^, while we are concerned with Eu¬ 
clidean processes, but we will see that nonetheless there is a close analogy. 
In both cases the semiclassical approximation of the Lipatov technique fails 
to give the correct dependence on (large) external momenta. The point in 
our case is that the coupling constant g of the theories of interest is dimen¬ 
sionful {g ~ {massy which has the consequence that skeleton graphs 
for processes with large external momenta p behave like {g'^/p)^ in Nth or¬ 
der perturbation theory, multiplied by a factorially-divergent factor. Just as 
in the Ringwald-Espinosa problem, the semiclassical level of Lipatov (even 
in the massive theory, which at least has a soliton) cannot reproduce the 
correct dependence on external momenta, which would require quantum cor¬ 
rections at least to Nth order to the semiclassical approach. That is, unless 
all the external momenta vanish|Q, the semiclassical Lipatov approach can 
only yield dependence on the external momenta which is a product of factors, 

^This subject, and the closely-connected subject of factorial divergences in perturbation 
theory, engendered a great number of papers and much controversy a few years ago. The 
upshot is that gauge-theory scattering processes in d=4 remain weak even at perturbative 
orders N :$> but that at high temperatures such theories can become strongly- 

coupled. For the author’s view, see Ref. 0; for recent summaries, see Refs. 
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one for each momentum, which is generally incorrect. In any case, a super- 
renormalizable theory with only one dimensionful parameter has no small 
dimensionless parameter, and is necessarily strongly coupled, invalidating a 
semiclassical approach. 

It is clear that d=3 gauge theories are strongly-coupled at small momen¬ 
tum, because they have no perturbative mass and /p grows large. We 
know, of course, that these theories generate 1,0 non-perturbative masses 
M ^ and the small-momentum expansion parameter is g'^/M = 0(1). It 
might be thought that the question of large-momentum perturbation theory 
is uninteresting because the Nth order corrections vanish rapidly. Unfortu¬ 
nately, there are 0(iV!) graphs, and so perturbative terms become large when 
N > p/g^- Another confounding factor for d=3 gauge theories is the lack of 
Borel summability, which makes it uncertain how the large-problem is to 
be solved by some sort of unitarizationj^ or other technique going beyond 
perturbation theory. 

We would like to study these questions for the gauge theories mentioned 
above, but this is too hard to do (at least for now), so we turn to a stand-in 
theory: This theory is, of course, irretrievably sick because its Hamil¬ 

tonian is unbounded below, but that will not bother us in the present in¬ 
vestigation. It is related to asymptotically-free theory in somewhat the 
same way the d=3 gauge theories are related to their asymptotically-free d=4 
counterparts; 05 is infrared-unstable and its coupling g also has dimension of 
(mass)^A We therefore expect to hud the same phenomena discussed above 
for gauge theories, including the lack of Borel summability because of Borel 
poles on the positive real axis. 

We attack the problem with direct graphical summation techniques^], 
which have proven useful for similar problems in theories with dimensionless 
couplings, such as 0| and 0g. Specihcally we will study, in 05 , the three- 
point one-particle irreducible skeleton vertex with bare internal propagators 
and vertices, or just the vertex, for short. In any held theory, the basic 
graphical building blocks are the skeleton graphs for vertices and propaga¬ 
tors, although they are little-studied at large N. These graphs have several 
special features, in comparison to the well-studied [jT^ simple ladder graphs. 
As functions of the external momenta they behave as naive power-counting 
suggests, in the massless theory; unlike, e.g., ladder or planar graphs, their 
multiplicity is roughly A^!; at large N the properties we will need are sta¬ 
tistically distributed with well-dehned averages and small huctuations; and 
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finally the bare skeleton graphs in theories with dimensionless couplings are 
the ones that correspond precisely to the semiclassical Lipatov approach^], 
even at asymptotically-large momentum. 

It is this last feature which is not true (at least in any way obvious to 
the author) in theories with dimensionful couplings, such as studied here. 
While there is certainly a connection to the Lipatov approach in the mas¬ 
sive theory, our graphical techniques correspond not just to the semiclassical 
hrst approximation, but instead to high-order corrections which instate the 
correct dependence on external momenta. 

The connection we hnd to the Lipatov soliton is a subtle one, and is 
best illustrated by looking at the results we will derive below for the Nth 
order contribution Ltv to the vertex, with the square of all three external 
momenta given by the single quantity p^\ 

Tn = N\{ag^)^{M^ + ( 1 ) 

(In this and subsequent equations the = symbol, to be used only at large iV, 
means that we ignore non-leading terms, constant factors, and factors of the 
type where 6 is a hxed constant.) Both a and a in (1) are positive numer¬ 
ical constants which we will estimate. When but p^ -C NM"^ it ap¬ 

pears that the large-p momentum dependence is of the form exp(—ap^/2M^), 
but this should not be construed as characterizing the behavior of the sum 
over N of all skeleton graphs, which also involves N such that iV 3> 1 and 
p^ 3> NM‘^. Except for this factor, when N becomes large at hxed (large) 
p, the momentum dependence goes away. Roughly speaking, the reason for 
this fading away of external momentum dependence at large N in the mas¬ 
sive case is that the graphs have 0{N) lines through which the momentum 
“diffuses”, but every one of these lines has the mass in it. The result is that 
when N is so large that the momentum term in (1) can be dropped, the sum 
of bare skeleton graphs beyond this value of N is correctly predicted by the 
Lipatov technique used at zero external momentum. From (1), the terms of 
this sum are N\{ag^/M)^, and one reads off the action of the Lipatov soliton 
as M/ag’^. As we have said before, there is no such soliton (in either (f)\ or 
d=3 YM) in the massless theory, where from (1) the behavior is essentially 
(3iV/2)!(pVp)^. 

There are two ways of estimating the action of the Lipatov soliton; One is 
direct numerical calculation (there are no analytic solutions), which we give 
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below; the other is to estimate the constant a from our statistical approach 
to skeleton graphs. This second way disagrees with the hrst by about 40%, 
which is comparable to the error in the graph-statistical approach as used 
earlier for theories with dimensionless couplings |^. The graphical estimates 
of the constants a, a require a non-trivial extension of this previous work. It 
involves graph-theoretic combinatorics which the author has been unable to 
hnd in the graph theory literature (as well as some known combinatorics, for 
which we give simple physicists’ arguments). 

Clearly, the next question is how to sum over iV in (1). A formal Borel 
integral can be given for the sum, by dropping the A^! in (1) and doing the 
resultant sum to hnd the Borel kernel. We will show how to express (in the 
= sense of (1)) this sum as an integral, and the formal Borel sum (which is 
the full skeleton vertex less the bare vertex) is a double integral of the form: 

_ poo ^ p 

V T^ = X / dudtt{—Y^^e-^-^^y+^'> [1 - (2eM)i/%-“xt]-h (2) 

Jo u 

In equation (2), x = ag^/M and y = ap^/M'^] the Borel-transform variable 
of integration is t. We will straightforwardly show that the singularity in 
t associated with infrared instability occurs at t = /M, the inverse of 

the classical Lipatov soliton, independently of the external momentum. This 
is, of course, what the Lipatov technique would give in the semiclassical 
approximation; what we get in addition by direct graphical techniques is 
a picture of how the summed vertex graphs depend on external momenta. 
That the Borel singularity is associated with classical solitons in d=4 gauge 
theories is an old argument of’t Hooft|^, which we now see is also true for 
lower-dimension infrared-unstable gauge theories. 

Two more questions now arise: (a) How are we to interpret these results 
for massive as having to do with d=3 gauge theories, which are massless in 
perturbation theory; (b) How is the Borel singularity to be regulated? Both 
answers must necessarily be speculative, without independent investigation 
of the gauge theories; we suggest some plausible but unconhrmed answers. 
As for the question of mass, we know that the d=3 gauge theories do generate 
a gluon mass M ~ non-perturbatively. We can roughly account for this by 
replacing the free propagators of the bare skeleton graphs by fully-dressed 
propagators, using a spectral representation. The result is that when the 
dressed propagators become massive, the quantity M in (1,2) is replaced 
by a quantity of order unity times the “physical” gluon mass. So we guess 
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that d=3 gauge theories do have graphical expansions, corrected for mass 
generation, somewhat like (1,2); of course, the “expansion” parameter g‘^/M 
is just a pure number which can only be calculated non-perturbatively. 

As for regulating Borel singularities, we speculate that this is done much 
the same way as is known to happen|^ in other models, where the semiclas- 
sical Lipatov solution (actually a sphaleron-like solution for gauge theories) 
conspires with the Borel sum of perturbation theory to remove unphysical 
imaginary parts. This amounts to replacing in (2) by e“^ —where 
t^iu^x) is the position of the pole in t in equation (2). We will discuss this 
further in Sec. 5. The result is that the regulated sum corresponding to (2) 
is strongly-coupled because ag'^/M is 0(1), but strong coupling has nothing 
to do with the external momenta. This is because the mass regulates the in¬ 
frared singularities of the massless theory, and the Lipatov soliton regulates 
the problems mentioned above which could occur when g^ but p < Ng‘^. 

In the regulated sum corresponding to (2) it is straightforward to show 
that the leading asymptotic behavior at large p is ~ just as one would 
expect from one-loop corrections and the naive argument that higher-order 
graphs are suppressed by more powers of p“^. 

Although there is no obvious way to associate (1) or (2) directly with a 
semiclassical Lipatov analysis, one gets the feeling that there really is some¬ 
thing semiclassical behind the direct graphical approach. This is certainly 
SO0 for theories with dimensionless couplings, where the skeleton graphs 
so far analyzed are essentially independent of momentum even when these 
are large, and therefore are in effect zero-momentum graphs, where the Li¬ 
patov technique surely works. And even for our theories with dimensionful 
couplings, the graphical analysis with its emphasis on well-dehned average 
properties of a very large number of graphs, and small fluctuations around 
these averages, suggests some underlying mean-held theory for their descrip¬ 
tions. This theory is presently unknown, and it would be very helpful to hnd 
it. 

2 Nth Order Bare Skeleton Vertex in (l)\ 

The Euclidean action is: 

S = j d^x\^di(t){xf + M‘^(t){xf + ^0(x)^]. (3) 


6 


We will be concerned for the most part with the behavior of the bare skeleton 
vertex function, which is dimensionless when an overall power of g is removed. 
The vertex at Nth order of perturbation theory, which we denote Tjv, involves 
divided by the Nth power of some quantity with the dimensions of 
mass. This quantity will certainly be the mass (if there is one) when all 
external momenta vanish, but it may involve external momenta too, when 
these are large compared to M. Exactly what happens depends on the graphs 
considered. Vertex graphs which are two-particle-reducible, such as the graph 
of Fig. 1, in general depend on both. For example, let pi vanish, and 
P 2 = P 3 = 3> M^. One simply approximates the propagator joining 

the vertices carrying momentum p by 1 /p^ to hnd that this graph is the 
product of this factor and a self-energy graph at zero momentum, which 
can only depend on g'^/M. More elaborate versions of this sort of behavior 
were studied long ago[PT![]. All of these older studies of graphs emphasize 
infrared singularities arising from highly-reducible graphs, such as sums of 
ladder graphs, but there is only one such (bare) ladder graph at every order 
N, so we cannot hope to study factorial behavior there. 

Bare skeleton graphs, however, have factorial multiplicity, and should 
correspond to the perturbative terms found by the Lipatov analysis. But 
at the semiclassical level this analysis decouples almost completely from the 
external momenta, as we will discuss in Sec. 4. To hnd the momentum 
dependence even approximately requires a graphical analysis. 

Consider an Nth-order vertex skeleton graph with bare vertices and prop¬ 
agators, an example of which is shown in Fig. 2. Such a graph, in graph- 
theoretic language, has girth four, which means that no loop in it has less 
than four lines; it has N loops, 3N lines, and 2N 1 vertices. It is both one- 
and two-particle-irreducible (except trivially by cutting the lines incident on 
any of the three external vertices). In any graph has the value, expressed 
in terms of Feynman-parameter integrals after the momentum integrals have 
been done, 

r(G) = S(G)r(N/2)lN^f J[dxU^u-^/^^/U + (4) 

Here 

[dx]3N = dxi... dx3Nd{l -^Xi), (5) 

U is the determinant of the graph, a sum of positive Nth-order monomi¬ 
als in the Feynman parameters Xj constructed according to rules reviewed 
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briefly below, and 0 is a sum of positive monomials of order + 1 in the 
Xj, multiplied by one or another of the pf according to rules reviewed below. 
The Mf are the squared masses of the internal lines; for later applications 
we distinguish them, but in the bare skeleton graphs each mass is just the 
mass M. Finally, S{G) is the symmetry number of the graph; at large iV, 
almost all (in the = sense) graphs have symmetry number of unity, and we 
will ignore this factor in what follows. 

The non-trivial difference between what we do here and previous work 
[P] is the appearance of the (p/U term in (4) to a power other than zero. 
For a dimensionless skeleton vertex in a theory with dimensionless couplings, 
this term either appears to the zeroth power or as a single logarithm (which 
can be eliminated by taking the Minkowski-space imaginary part). In such 
cases there is no dependence on external momenta at all, and these might as 
well be zero. This means that the sum of (the imaginary part of) skeleton 
graphs at large order must reproduce the semiclassical Lipatov value, which 
is only completely correct when all the momenta vanish. But in the present 
case the p/U term is present and must be dealt with. This is, by the way, 
exactly what must be done to understand the perturbation-theory analog 
of the Ringwald-Espinosa problem, which in part stems from the fact that 
graphs with a large (0(1 /A)) number of external legs carrying momentum, in 
a theory with dimensionless coupling A, have an equally-large negative power 
of some combination of the external momenta. This negative power comes 
from the cp/U term. The semiclassical Lipatov approach does not give this 
dependence correctly. 

We will deal with this p/U term just as was done earlierP| for the case 
when it did not appear: Just replace p as well as U in (4) by their average 
value (averaged over the Feynman-parameter integral). Such an approach 
certainly cannot work for many types of Nth-order graphs, for example, the 
vertex ladder graph. Ladder and other highly-correlated graphical structures 
have (/-functions which can vanish very rapidly compared to U, and thereby 
generate power-law infrared singularities (that is, inverse powers of external 
momenta which would naively be expected when the momenta are large 
get replaced by inverse powers of the mass M; this is how Regge poles are 
generated |]T3 in p^). But large-iV skeleton graphs are very random, and one 
cannot make p vanish without making 0{N) Feynman parameters smallQ 

^When all the parameters of J loops are made small, U vanishes like the Jth power of 



the region of parameter space where this can happen is too small to lead 
to an infrared singularity. This statement is just another way of stating 


Weinberg’s|^ theorem. So we feel justihed to replace both 0 and U by their 


parameter-space averages, as was earlier justihed simply for U aloneQ. 


The rules for U and 0 are well-known ||15||. To hnd U, hnd the chord sets, 


which are the sets of all lines which, if cut, turn the graph into one connected 
tree, called a spanning tree.. The monomials of U are the products of x* 
constituting the chord sets. To hnd 0, hnd the cut sets, or lines which, if 
cut, turn the graph into exactly two connected trees. Write the product of Xi 
for each cut set, multiply the result by the square of the sum of the external 
momenta howing into or out of either tree, and sum; the result is 0. 

As discussed above, we argue that for hxed and large order N the sum 
of all bare skeleton graphs of that order is (in the = sense) reasonably well- 
estimated by replacing U and 0 in equation (4) by their averages over the 
Feynman-parameter integral, and multiplying (4) by Qn, the number of 
graphs of Nth order. It is easy to compute the parameter-integral aver¬ 
age of any monomial; the two of concern are 0{N) and 0{N + 1) and we 
have 


r(3iv) 



( 6 ) 


T(3N + K) 


SO 



(7) 


(xi . . .Xiv+l) = • --Xn)- 


( 8 ) 


It follows that {U) is found by multiplying the result in (7) by the com¬ 
plexity Cat of the graph, that is, by the number of spanning trees. This 
number has been found [^j numerically earlier: 



(9) 


Then 




( 10 ) 


the smallness parameter, so there are cancellations in the ratio 0/C/. 
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The estimation of 0 is a little more complicated. First, we need the 
number of cut sets of the graph, and in the Appendix we show that this 
number is just 20^- However, this is not exactly what we need. Each cut set 
monomial is multiplied by the square of the momentum entering either of the 
two trees; sometimes this is one of the external momenta pj and sometimes 
it is zero (see Fig. 3). In the Appendix we give a sketchy argument that the 
fraction (3 of cut sets containing an external momentum is 1/2. Since each 
monomial in 0 has the same average, it is clear that (0) is the same if we 
replace each pj by the average squared momentum p^\ 


v^^^ivl + vl+vD- 

Then using these remarks and equation (8) we find: 

W = 2;3A(f/> = A(C/> 


( 11 ) 


( 12 ) 


Finally we need the elementary integral 



(13) 


and the number Qat of A^fh-order 0^ graphs. This number is known |]T6|, and 
we give a physicists’ derivation of it in the Appendix: 

Qk^IVkI)" (14) 

This number includes all graphs, skeleton or otherwise, but as the derivation 
in the Appendix shows, almost all the graphs are skeleton graphs. 

We now replace 0 and U in (4) by their averages, and multiply by Qa? 
to find an estimate for the sum of Nth-oidei skeleton graphs as given in 
equation (1), with 


217/2 

37(47r)5/2 


2.96 X 10■^ 



(15) 


Three remarks: First, at pi = 0 the graphical expansion should yield the 
Lipatov results. We study this in Sec. 4, finding for the coefficient a the 
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value 4.55 x 10““^. The graphical value is about 2/3 the Lipatov value, and is 
less than it, in line with general arguments 0 based on Holder inequalities. 
Second, one might compare the above value for a with that found by graphical 
techniques for 0g, where a is dehned by the expansion 

T^ = N\{agy (16) 

and the dimensionless coupling qq appears in the action just as it does in 
equation (3). Dehne the dimension as d = 6 — e; then the present hve- 
dimensional results and the earlier P| six-dimensional results can be written 
in one formula as: 

2ll-5e/2 

“= 38 -./ 2 ( 4 ^) 3 -./ 2 - ( 1 '^) 

This is a kind of e-expansion which can be used for infrared-unstable theories. 
Third, when one goes beyond bare skeleton graphs by dressing the vertices 
and propagators in such a graph, one will hnd logarithmic terms which pre¬ 
sumably will change powers of momenta and lead to anomalous dimensions. 
This is a subject for future analysis. 


3 The Formal Borel Sum 

We have now established the result (1) (repeated below for convenience) 

Tjv = + ap^N)-^/^ (18) 

and would like to know how to write it in the usual Borel-integral form. Of 
course, this integral will have a singularity on the positive real axis, because 
all the terms of (18) are positive. We will discuss regulation of this singularity 
later. 

Given a formal sum 

h{x) (19) 

the Borel prescription is 

poo 

h{x) = / dte~^g{xt) (20) 

Jo 

where 

( 21 ) 
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So from (18) we need the sum 

9{x, y) = Y^ a;^(l + y/N)-^/^; x = ag^/M, y = ap^/M^. (22) 

1 

This non-standard sum can be written (in the = sense) as an integral, using 
the identity 

+ iv)-^/2 = (23) 

Then Stirling’s formula gives 

(1 + y/N)-^/"^ = (2e)^/2 (24) 

Use this to do the sum in (22): 

poo 

g{x,y) = x due~'^'^y ~^^\—)^'^^[1 — a;(2eM)^/^e“'“], (25) 

JO u 

and from (20) the formal Borel sum of (18) is equation (2), repeated for 
convenience: 

poo Op 

J2 Tn = xJ^ dudtt{—y/^e-*-^^y+^'^[l-{2euy/‘^e-^xt]-\ (26) 

The singularities of this double integral are endpoint branch cuts, express¬ 
ing a dependence on p rather than and a Borel pole where the denominator 
in (26) has a second-order zero, or pinch singularity. (A single zero not at 
an endpoint can be evaded by shifting the contour.) The second-order zero 
occurs where the function (2eM)^/^e““ has a maximum, which happens at 
u = 1/2, and the maximum value is unity. Therefore, the Borel singularity 
is at to = = M/ag'^, and the Borel ambiguity is e~^° = . This 

is, as expected, the exponential of a classical soliton action associated with 
the Lipatov technique. So we should be able to hnd the same action (he., 
the same value of a) by looking for classical solitons of (p^. This we do in the 
next section. 
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4 Lipatov Analysis and the Need for Mass 

We briefly review the Lipatov technique, showing how it fails at the semi- 
classical level to get the right dependence on external momenta. Then we 
show that the technique fails, even at zero momentum, for zero mass, and 
hnally we hnd the appropriate soliton for non-zero mass. 

Write the Nth order contribution to the full (reducible) vertex in path- 
integral form: 

Z~^ jJ J d^x exp{ipi ■ Xi)(j){xi) (27) 

i 

where Sq is the free action. At large N the external wave-function factors are 
to be evaluated by inserting the classical solution 4 >c{xi — a) and integrating 
over the translational collective coordinate a; the result is a product of three 
factors each depending on only one external momentum. This is clearly 
wrong for the skeleton graph of Fig. 2, unless all the momenta vanish. But 
if all the momenta vanish a mass is necessary to construct a dimensionless 
expansion parameter g'^/M. The smaller the mass the more strongly coupled 
the theory is, but in an infrared-unstable theory massless at the Lagrangian 
level (including 01) there is a lower bound on the mass, arising from tachyonic 
poles in the zero-mass limit of such theories |0, |17| . 

The Lipatov technique understands this non-existence of the massless 
theory at zero momentum as the non-existence of a massless soliton with 
hnite action (of course, the massless soliton is precisely what does arise in 
theories with dimensionless couplings). The classical equations coming from 
the saddlepoint of the action^ in equation (27) above is: 

□0 - M2<F + = 0; K = [ d^xcf)^. (28) 

K J 

Let us look for a spherically-symmetric solution for the massless case. Even 
though the underlying held theory has a dimensionful coupling constant, 
it is absent in (28) and this equation is independent of the scale of length 
chosen (there is a collective coordinate for scale). Call this scale p and change 
variables to: 

^ = e-h (29) 

^We can now drop the external wave-function factors. 
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(Note that from (28) K ~ so 0 has its canonical mass dimension of 3/2.) 
Performing this change of variables in (28) yields 

+ -^(1 = 0- (30) 

The dots indicate derivatives with respect to t, so this is the equation of 
motion of a particle of coordinate q with friction in a potential 

q^ ( 7 ^ 

= ( 31 ) 

This potential has a relative maximum at g = 0 and a minimum at g = 2; it 
reaches the value zero at g = 0, 3 so that without friction there would be a 
soliton going between these two points. But this is impossible with friction, 
and any solution vanishing at x = oo {t = —oo) necessarily settles to the 
minimum at g = 2. Indeed, g = 2 or from (29) 0 ~ is an exact solution, 
but with inhnite action. 

There is an analog of this in d=3 YM theory. The massless version 
has only one spherically-symmetric soliton, the Wu-Yang soliton, and it has 
inhnite action too. Indeed, one can repeat the steps of writing down the 
equations and changing variables rather as in (28, 29) to hnd the same sort 
of motion-with-friction analog. So the Lipatov analysis fails here, too, as we 
know it must. 

Now turn to the massive case, and write 

KM"^ 

m = (32) 

The equation for p{u) is: 

p" + ~ + T(pi - P) = 0 (33) 

and it does have a hnite-action solution (with 0 approaching a constant at a; = 
0 and vanishing exponentially at inhnity). We have determined the solution 
and the action numerically. One can then easily evaluate the coordinate 
integrals in (27) and thereby hnd the semiclassical Lipatov-technique value 
for the contribution Tjv to the Nth-order vertex at zero momentum: 

r^ = Y!(^)'^; a = 4.55x10-". (34) 
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This value of a is to be compared to the graph-technique value of a = 2.96 x 
10“^ given in (15). 

Again, the existence of a massive soliton has its analog in d=3 YM theory, 
which also has one||18||. Its action is dvr x bAM/g'^. Of course, the mass 


does not appear in the clasical action, but is generated by non-perturbative 
quantum effects^, 0 which can be described in terms of an effective massive 
action, which does have a soliton. Even thought M ~ the mass vanished 
to all orders of perturbation theory. So although it is essential to include 
the mass in order to stabilize the theory, once one does so in, e.g., the bare 
skeleton graphs one has really gone beyond perturbation theory altogether. 
In the next section we speculate on this and other non-perturbative features. 


5 Beyond Perturbation Theory 

Aside from the question of how to incorporate a non-perturbative mass in 
skeleton-graph expansions like (1) or (2), there is the problem of regulating 
the Borel singularity in (2). 

In 05, if not in d=3 gauge theory, it is straightforward to incorporate 
dressed propagators into the bare skeleton graphs. One writes for the prop¬ 
agator on line i a Lehmann representation; 

A(Pi) = [ / dMfp(A4) = 1 (35) 

Pi ' J 

SO that the right-hand side of equation (4), which gives the contribution of 
the iVth-order graph G, becomes 

(36) 

As before we will replace functions of the Feynman Xi by their averages, and 
(xj) = l/3iV, so the J2xiMf term in the denominator is replaced by the 
average value cr: 

(V) 

To the extent that the positive weight functions p{M‘f) behave like probability 
distributions we can invoke the central limit theorem to replace a by the 
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average of Mf over p, clearly a quantity differing from the physical mass 
by a factor of order unity. So our main results (1,2) continue to hold, with 
dehned as this average. 

It is somewhat delicate to apply this line of argument to d=3 gauge theory, 
for several reasons. First, one cannot maintain gauge invariance without 
dressing vertices as well as propagators; second, even if the internal vertices 
are dressed in the skeleton graphs, the hnal result is gauge-dependent and 
must receive corrections such as used in the pinch technique[|^, |T^; third, 
in any theory where there is no mass at the classical level, the quantum 
mass must depend on the propagator momentum and vanish at large p^. 
It is known poll that this running mass vanishes like g‘^{G‘fj)p~‘^. However, 


presumably one should not replace in (1) or (2) by M‘^{p‘^), but rather 
by M‘^{p‘^/N), just as the graphical analysis tells us to do for the kinetic- 
energy terms of propagators. 

Arguments like this aside, there is another reason to believe that it may 
be a decent approximation for d=3 gauge theory to use (1,2) with the zero- 
momentum quantum mass M in it, thus going beyond perturbation theory. 
This is that the Lipatov-technique result at zero momentum is likely to be 
the right starting place to understand large-order dressed skeleton graphs, 
and this technique will incorporate whatever quantum mass is generated in 
the full theory. 

Now turn to the problem of regulating the Borel singularity. It is known 11^ 
in various quantum-mechanical and held-theoretic problems that the Borel 
singularity, which may give rise to an imaginary part having no physical in¬ 
terpretation, is cancelled by the contribution of solitons such as the one we 
have found. We do not attempt to study this problem in detail here, but 
simply suggest a plausible answer, which is to subtract off an exponential 
term in the Borel integral which gets rid of the singularity. So (2) or (26) 
gets modihed to 

_ poo Op 

yr^= / dndtet(—)^/2g-n(y+l)r (38) 

Jo u 


where 


to(a:,n) = 


(39) 


x{2euYI‘^' 

It can now be easily checked that there is no sign of a large contribution from 
large momentum p, as one might have naively expected from a behavior like 
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N\{g‘^/p)^ when N is large and p ~ Ng"^. The leading large-p behavior is 
what is expected from one-loop corrections, namely g‘^/p. 


6 Conclusions 

In this paper we have studied the large-order perturbative behavior in theo¬ 
ries where external momenta are expected to play a signihcant role. The two 
cases of interest are in theories where the coupling constant is dimensionful 
(such as d=3 gauge theory or or, in d=4 gauge theories, the case of multi¬ 
particle amplitudes where the number of particles is 0(l/a)@, ||, The 

troublesome point is that the widely-used Lipatov technique cannot sense 
these external momenta in the semiclassical approximation. To see their ef¬ 
fect one has to do something much like the graphical analysis we give here, 
or, much better, hnd a Lipatov technique based on some sort of effective held 
theory which allows incorporation of the external momenta. 

Even at zero momentum in d=3 gauge theory, one must use solitons of 
an effective action incorporating a non-perturbative gluon mass. Otherwise 
the Lipatov technique falls apart, for lack of a massless soliton with hnite 
action. Strictly speaking perturbation theory at large order never makes 
sense in this theory, because we have found that both small-p and large-p 
problems (the latter coming from large N) are resolved by a hnite mass, 
which is absent in perturbation theory. In massless theory at large p 
one can always go to N so large (of order {p/g‘^y^^, based on (1)) that 
perturbation theory breaks down, quite the opposite of what one expects 
in a superrenormalizable theory, where higher orders are supposed to fall 
oh more rapidly in p. However, in the massive theory as N increases the 
momentum dependence is overtaken by a mass dependence such that the 
expansion parameter becomes g'^/M. Perturbation theory still breaks down, 
but because of the infrared singularities of the theory. 

It would, of course, be valuable to extend our graphical work to such a 
gauge theory, but this seems to be very hard indeed at present. Again, an 
ehective-action approach would be much better. 

The graphical analysis goes beyond earlier large-work in including the 
ehect of dimensionful denominators through the 0 function of Sec. 2. It 
would be interesting to extend this work to problem of multiparticle am¬ 
plitudes in d=4 renormalizable theories, to see if the kind of random-graph 
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arguments we use here and earlier [Q can be made to work. What we have 
done is a minimum hrst step. In view of the mass dimensionality of d=4 
scattering amplitudes, which is 4 — iV for an Wpoint vertex, a behavior like 
(A/p)^ in a theory with dimensionless coupling A corresponds to only that 
part of the high-energy hxed-angle amplitude in which the order of pertur¬ 
bation theory is correlated with the number of legs. More work is needed to 
study the general case. 
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Appendix: Graphical Combinatorics 


Here we discuss three questions: (a) The number of Nth-oidei 0^ graphs; 
(b) The number of cut sets at this order; (c) The fraction of cut sets which 
carry hnite momentum. 

To count graphs, we can construct an A^ + Isf-order graph by adding one 
line to any A^fh-order graph in all possible ways. There are (3A^ + 3)^/2 
ways to do this (we never consider graphs with one-line loops, and the added 
line can begin or end on any of the 3N internal lines or the three external 
lines). But this overcounts by a factor of 3N -|-3. To see this, mark any other 
line in the Nth-order graph (not the one added to make the N -\- 1st order). 
Exchange the marked line and the added line. This gives the same graph in 
N -|- 1st order, but constructed from a different graph at Nth order. So the 
relation between the number of graphs Qn+i at N 1st order to Qat is 


Qn+1 — 2 ^^ 


On = 


(40) 


Almost all of these graphs are skeleton graphs. Let the A^fh-order graph 
be a skeleton. Then, having chosen one of 3N -|- 3 positions to place one 
end of the added line, there are only 3N ways to place the other (it cannot 
go on the line it started, or on lines separated from it by only one vertex). 
Continuing as above, one hnds that Qn for skeleton graphs differs from that 
for all graphs by a factor N'^, where c is a hxed constant. So in the = sense 
the two counts are the same. For a full discussion of the graph-counting 
problem, see ref. [jK . 

Next we count the cut sets, which enter the dehnition of cp. A good 
heuristic for this is to start from the complexity, or number of spanning 
trees. The complexity Cat is easily bounded: To form a tree, one must cut 
N out of 3N lines, so 


Cn< 


{3N)\ 

{2N)\N\ 



(41) 


. The actual number, computed numerically0, is (16/3)-^. The Nth roots 
differ by only 25%. The true number is smaller than 27/4 because not every 
set of N cuts yields a connected graph with no loops. 

We can End the cut sets by cutting the tree one more time; that is, a cut 
set comes from the original graph by making A^ -|- 1 distinct cuts. Thus the 
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number of cut sets Cat is bounded by: 


Cm< 


(3iV)! 

(iV + l)!(2iV- 1)! 




This suggests that Cn = ‘^Cn, for which we now argue in the same way we 
found Qi\f. Take a spanning tree and mark any line not already cut; there 
are 2N ways to do this. Now exchange the marked line with any line already 
cut. Any of the N ways of doing this exchange repeats a graph, just as we 
argued for Qat. Therefore we divide 2 N by N to come to Cat = 2Cn- 

Not in all cases does one of the two trees have external momentum flowing 
into (or out of) the tree (recall Fig. 3). We estimate this fraction from a 
particular set of graphs of the type in Fig. 4, in which N vertices he on 
the top line and another N on the bottom; these are joined in all possible 
ways by lines which we call crossing lines. There are clearly A^! such graphs, 
less than Qm, but we have no reason to believe our argument fails for more 
complicated graphs. 

First we make a spanning tree by cutting the line 1-2 a certain number of 
times, the line 1-3 a certain number of times, and the crossing lines a certain 
number of times; an example is shown in Fig. 5, in which two of the crossing 
lines remain uncut, and there are J cut crossing lines between vertices 1 and 
A, as well as K cut crossing lines between B and 3. We need to cut one more 
line to make a cut set. Evidently, if we cut any of the J (noncrossing) lines 
between vertices 1 and A, or any of the K lines between B and 3, we form 
a graph of two connected trees both of which have a total of zero external 
momentum flowing into the graph; we call this a null cut. Any other cut 
leads to trees carrying one of the external momenta pi] we call this a non¬ 
null cut. Elementary combinatorics which is not worth recording here shows 
that, on the average when N, J, and K are all large, the ratio of null cuts 
to non-null cuts is just L/(L-|-2), where L is the number of uncut crossing 
lines. This ratio approaches unity as L becomes large, so we conclude that 
about half the cuts are non-null and enter into the expression for 0. To 
form (j) all the Feynman parameters in a non-null cut are multiplied together 
and then multiplied by the square of the external momentum flowing into 
either tree. As explained in the text we replace 0 by its average wherever it 
appears. Since the average of every monomial of A^ -|- 1 Feynman parameters 
has the same value, it is evident that (0) is a symmetric function of the pj, 
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and can be calculated by multiplying the monomial average by the number 
of non-null cut sets and by the average squared momentum 


P 


2 


1 

3 


{pI + pI + pD- 


(43) 
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Figure Captions 


Fig. 1. A two-particle-reducible graph. 

Fig. 2. A bare skeleton graph; No internal propagator or vertex parts. 

Fig. 3. A graph cut into two connected trees. Neither tree carries net 
external momentum flowing into the graph. 

Fig. 4. A skeleton graph drawn by connecting N vertices on line 1-2 with N 
vertices on line 1-3. 

Fig. 5. A spanning tree for a graph of the type in Fig. 4. 
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